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Coherent oscillations of a scalar field can mimic the behavior of a perfect fluid with an equation-of- 
state parameter determined by the properties of the potential, possibly driving accelerated expansion 
in the early Universe (inflation) and/or in the Universe today (dark energy) or behaving as dark 
matter. We consider the growth of inhomogeneities in such a field, mapping the problem to that of 
two coupled anharmonic oscillators. We provide a simple physical argument that oscillating fields 
with a negative equation-of-state parameter possess a large-scale dynamical instability to growth of 
inhomogeneities. This instability renders these models unsuitable for explaining cosmic acceleration. 
We then consider the gravitational instability of oscillating fields in potentials that are close to, but 
not precisely, harmonic. We use these results to show that if axions make up the dark matter, then 
the small-scale cutoff in the matter power spectrum is around 10~^^ M^. 



I. INTRODUCTION 

Observational evidence for inflation, for accelerated ex- 
pansion in the current Universe Jj, and for nonbaryonic 
dark matter has motivated the search for exotic forms of 
energy. Scalar fields have been thoroughly investigated 
for both pressureless and negative-pressure matter. If the 
kinetic and potential energies of the scalar field are nearly 
equal, then the field behaves as pressureless matter. The 
most important example is axion dark matter, which can 
be described as coherent oscillations of the axion field 
in a nearly harmonic potential. Alternatively, if the ra- 
tio of potential energy to kinetic energy is sufficiently 
large, cosmic acceleration can be induced. Quintessence 
achieves this with a rolling scalar fieldjthe rolling slowed 
by Hubble friction • Spintessence [1, |3] achieves this 
with a complex scalar field rotating in an internal C/(l) 
symmetric potential; here, the centripetal acceleration, 
rather than Hubble friction, prevents the scalar field from 
falling directly to its minimum. A coherently oscillating 
scalar field can also drive cosmic acceleration. 

Coherent oscillations in a harmonic potential behave 
as nonrelativistic matter, but oscillations in a more gen- 
eral potential can mimic a perfect fluid with an arbitrary 
equation-of-state parameter. For example, in the pres- 
ence of a power-law potential V{(j)) oc |0|", one flnds that 
w — {n — 2)/{n + 2) The power-law index n de- 
termines, through the virial theorem, how the kinetic- 
and potential-energy densities are apportioned over one 
oscillation cycle. When n < 1, we have w < —1/3, 
illustrating that coherent oscillations of a scalar field 
might drive a period of inflation in the early Universe 
111 T3I or p rovide a candidate for dark 
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The purpose of this paper is to investigate the dynani- 
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ical (i.e., those arising from the scalar- field dynamics) 
and gravitational instabilities to the growth of inhomo- 
geneities that may arise in oscillating-field matter. The 
presence of such inhomogeneities are important for de- 
termining the viability of these models for describing in- 
flation, dark energy, and dark matter. 

The question of dynamical stability of oscillating-field 
dark energy has been analyzed, in the context of inflation, 
in Refs. [20, [Ml, concluding that oscillating fields that 
give rise to accelerated expansion are indeed unstable to 
the growth of inhomogeneities. This prior work analyzed 
the equations of motion, and it focused on the resonant 
growth of perturbations on small scales. Ref. [l^] con- 
sidered further the possible nonlinear evolution of the 
field. Refs. [1, [ij, [3| speculated that accelerating oscil- 
lating potentials may have a large-scale instability, but 
they did not carry out a full instability analysis. In 
spintessence, an analysis of the linearized equations of 
motion for scalar-field perturbations about the homoge- 
neous solution shows that many models driving cosmic 
acceleration (for a power-law potential, large-scale insta- 
bilities set in for n < 2) are dynamically unstable [1, [2^ . 

In addition, prior work [13, [IE] has shown that gravita- 
tional instabilities of oscillations in a harmonic potential, 
suitable for oscillating-field dark matter, are suppressed 
on sufficiently small scales. This scale determines the 
small-scale cutoff in the matter power spectrum, and has 
potentially interesting implications for cosmology. 

Our principle new contribution to the dynamical sta- 
bility of oscillating-field matter is a simple physical pic- 
ture of the criterion for stability at large scales. We show 
that the perturbed scalar-field equation of motion is iden- 
tical to that of two coupled anharmonic oscillators. This 
picture then allows us to provide a simple understand- 
ing of why oscillating potentials with negative pressure 
should be unstable at large scales, while those with pos- 
itive pressure should have large-scale stability. We then 
verify these conclusions analytically for potentials that 
are nearly harmonic (i.e., nearly pressureless) and numer- 
ically for a broader range of potentials. Our gravitational 
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analysis generalizes prior work [2J, |25| by considering po- 
tentials that are nearly, but not precisely, harmonic. 

Below, we first review (in Section II) the homogeneous 
evolution of an oscillating scalar field. We then provide 
in Section III a heuristic discussion of the dynamical in- 
stability, beginning first with an explanation for the ori- 
gin of the dynamical instability in spintessence models. 
Section IV shows analytically that for nearly harmonic 
potentials, instability occurs for negative-pressure poten- 
tials and stability occurs for positive-pressure potentials. 
Section V verifies this conclusion numerically for more 
general potentials. Section VI includes gravity in the 
analysis and works out the gravitational-instability scale 
for nearly harmonic potentials. We discuss here the ap- 
plication to axion dark-matter models, working out the 
small-scale cutoff in the matter power spectrum. Section 
VII presents some concluding remarks. 



II. HOMOGENEOUS EVOLUTION 

Consider a scalar field with potential V{4>) = V"(— 0) 
minimized at V{(j) = 0) = 0. The equation of motion for 
the scalar field is 



4> + v\(i)) = o, 



(1) 



where the dot denotes derivative with respect to time t 
and the prime denotes a derivative with respect to the 
scalar field 0. The scalar field will undergo periodic mo- 
tion in this potential, with a period 
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(2) 



that depends, most generally, on the scalar- field ampli- 
tude (p q. T he period can also be written in terms of the 
action [l^l 



J 



# v/2[F(</>o)-y (</>)], 



(3) 



as T = dJ/dVo, where Vq = V{(f>Q). 

For V{(j)) oc |0|", the angular frequency w(0o) ^ 



|(l/2)-(l/n) 



Thus, for n = 2, the frequency is 



amplitude-independent. For n > 2, the frequency in- 
creases with amplitude, and for n < 2, the frequency 
decreases with amplitude. 

On timescales long compared with the oscillation pe- 
riod, the scalar-field oscillations behave like a perfect 
fluid with energy density p and pressure p, 



p = {^0'/2 + V 

p = Uy2-v 



(4) 
(5) 



where the angle brackets denote a time average. The 
virial theorem tells us that l(jP/2\ — (1/2)(1 -I- w)Vo and 



(1/(0)) = (1/2) (1 - w)Vo, where Vq = V{<i)o) and w = 
pI p — (n — 2)/{n -I- 2) is the equation-of-state parameter. 
These results can be summarized in terms of the action 
by [il 
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(6) 



Differentiating the relation T = dJ/dVo and re-arranging 
algebraically, we obtain an expression for the variation of 
the frequency with field amplitude. 
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(7) 



Thus, the frequency increases with amplitude unless 
Vodw/dVo +w <Q. 

If w is independent of (^Oj as for power-law potentials, 
then the frequency will decrease with amplitude when 
w < and increase with amplitude when w > 0. If w is 
allowed to change with ^o, then we can have w < and 
duj/d(j)o < 0, but only temporarily. 

There is also a nice geometric interpretation of the 
value of the equation-of-state parameter Q. Using en- 
ergy conservation and the field equation, we can write 



V- 



(1 + 3w) 



Vo 



(8) 



The left-hand side is the average of the intercept of the 
tangent to the potential, and is in general positive for 
convex (about = 0) potentials and negative for non- 
convex potentials. A few simple example potentials are 
shown in Fig. [1] to illustrate this, although it is easy to 
imagine potentials with more complicated features. In 
general, to produce accelerated expansion, there must 
be a relatively flat region of the potential with positive 
energy somewhere along the oscillation cycle, according 
with the intuition that the energy density must be poten- 
tial dominated in order to produce accelerated expansion. 

We also note that the sound speed (squared) is in 
general different than the equation-of-state parameter w, 
since the sound speed is given by p^] 



, d(wp) dw 
ct = -^-^ = Vo — 
p dVo 



(9) 



We thus see that the sign of the sound speed is the same 
as the sign of dT/dVo- 



III. HEURISTIC DISCUSSION OF 
DYNAMICAL INSTABILITY 

A. Preview: the spintessence instability 

We begin by considering the growth of perturbations 
in spintessence. These models introduce a complex scalar 
field (f) with a [/(l)-symmetric potential V{(j)) = y(|0|). 
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FIG. 1: Examples of potentials. The red dashed potential pro- 
duce an equation of state corresponding to accelerated expan- 
sion {w < —1/3). The average intercept for some amplitude 
is shown as the red dot sd, V > 0. The harmonic potential, 
which produces an equation-of-state parameter ui = 0, is the 
blue dot-dashed curve, and its average intercept the blue dot 
at y < 0. Oscillations in the solid green potential produce an 
equation-of-state parameter w > 0, with an average intercept 
at the green dot a.t V < 0. 



tential, the angular momentum (per unit mass) increases 
as r increases. If the particle at large r runs ahead, as it 
will when y(|(/)|)/|(;i|-y"(|(/)|) < 0, then it pulls the inner 
particle forward, donating some of its angular momentum 
in the process. In this way, the inner particle evolves to 
a slightly larger r and the outer to a slightly lower r, 
decreasing their separation, and implying that the per- 
turbation is stable. If, however, the particle at larger r 
runs behind, as it will for V'{\(j)\)/\(l)\-V"{\(j)\) < 0, then 
it pulls on the inner particle, taking angular momentum 
from it. The inner particle must then evolve to a smallcr- 
r orbit and the outer particle to a larger-r orbit. In this 
way, the small initial separation between the particles is 
amplified, and the perturbations are unstable. Readers 
familiar with accretion-disk physics will recognize this in- 
stability as the source of angular-momentum transport in 
disks (see, e.g., Ref. ^) . 

For F(|(/)|) cx |0|", the instability sets in for n < 2. It 
is thus concluded that power-law spintessence potentials 
with negative pressure are subject to this instability and 
are therefore unsuitable as dark-energy candidates. 



The scalar field moves in a circular orbit in the poten- 
tial at some constant amplitude \(j)\. There is kinetic 
energy associated with the spinning and potential en- 
ergy associated with the displacement of the field from 
the minimum; the balance between the two is such that 
the equation-of-state parameter is w = [|(/)|y(|0|) — 
"^VilM/MV'M) + 2V{\(j)\)]. For example, if VM) cx 
|(^|", then w — {n — 2)/{n + 2), matching the result for 
an oscillating real field. 

Ref. 3 showed (see also Refs. [26|) that this coherently 
spinning field remains stable (neglecting gravity) to small 
perturbations if 1/'(|0|)/|(?!)| - F"(|(/>|) < and unstable, 
at sufficiently long wavelengths, if (|0|)/|(?!)| -y" (101) > 
0. 

This result can be understood simply. The dynamics 
of perturbations in the linear regime are identical to the 
evolution of two particles undergoing circular orbits in a 
two-dimensional circularly symmetric potential V{r); the 
gradient-energy density in the scalar field acts as a spring 
of force constant (where k is the Fourier wavenum- 
ber of the perturbation) connecting the two particles. 
The radii of the two orbits differ initially by only a tiny 
amount. In the absence of any coupling between the two 
particles, they will evolve independently, spinning at two 
slightly different angular frequencies and each staying at 
its original amplitude. If F'(|0|)/|(/>| - ^"(101) < 0, then 
the particle at larger r will run slightly ahead, and if 
V'M)/\(t>\ - V"{\(j)\) < 0, then it will run slightly be- 
hind. 

Now suppose there is a very strong spring that attaches 
the two particles. In this case, the two particles will be 
bound to each other, and there will be no growth of per- 
turbations. If there is a very weak spring, there will be 
energy transfer between the two particles on a timescale 
longer than the period of oscillation. In any confining po- 



B. The oscillating- field instability 

A similar argument can be applied to the growth of 
perturbations in oscillating-field matter. Adding a per- 
turbation S(j){x, t) to the homogeneous solution for the 
scalar field, the linear equation of motion for the pertur- 
bation obtained from Eq. ([T|) is 

5(l)-\l'^5(j) + V"{(j))5(l) = 0. (10) 

In linear theory, each Fourier mode of the density field 
54>-^ evolves independently and satisfies an equation (sup- 
pressing the k subscript), 

5(t> + k^5(l) + V"{4>)5(t) = 0. (11) 

This equation of motion is identical to that for the sep- 
aration between two particles connected by a spring of 
force constant fc^ moving in a potential 1^(0) with a sep- 
aration 54> = 4>2 — 4>i 4>i,4>2- 

In this picture, particle 1 is released at rest from some 
initial height ipo (where we will later identify 0o as 
the amplitude of the background oscillations), and par- 
ticle 2 is released slightly higher at (j)2 — ipi + ^4>- Be- 
cause the frequency of oscillation will in general be am- 
plitude dependent, if there is no spring connecting the 
two masses, then each oscillates at its own frequency. If 
^'{4>o) > 0, then the higher-amplitude mass (particle 2) 
runs ahead of the lower- amplitude mass (particle 1), and 
vice versa if cj'(0o) < 0- 

If the spring connecting the two masses is extremely 
strong, then the two masses oscillate together. However, 
if the spring is weak, then energy can be exchanged be- 
tween the two particles. If (jj'{(j>o) > 0, then particle 
2 runs ahead, pulls on particle 1, and consequently do- 
nates some of its energy. Particle 1 then moves to a 
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higher-amplitude orbit, particle 2 to a lower one, and the 
separation between them diminishes. Such perturbations 
are stable. If uj'{(f)o) < 0, then particle 1 runs ahead and 
donates energy to particle 2. Particle 1 thus moves to 
a lower-amplitude orbit, and particle 2 to a higher or- 
bit. In this case, the initial separation is amplified, and 
perturbations are unstable. 

We conclude that the sign of (^'{(po) provides a cri- 
terion for stability. This can be obtained from Eq. ^ 
as sign{oj') = sign{VQdw / (Wq + w). As discussed 
above, power-law potentials with negative pressure have 
^'(0o) < 0, and so will develop large-scale (small k) insta- 
bilities, rendering these models unsuitable for accounting 
for dark energy. Using the geometrical condition Eq. ([S]) 
as guidance, it is possible to find potentials that produce 
cosmological acceleration and do not develop large-scale 
instabilities, but only for a small range of amplitudes. 
The potential in the left cell of Fig. [51 which is of the 
form suggested by Ref. will exhibit long-term cos- 
mic acceleration, but only short-term stability as the am- 
plitude of oscillation inevitably decays. The potential in 
the right cell of Fig. [5] will exhibit long-term stability, 
but only short-term cosmic acceleration. The analysis of 
the viability of these models is somewhat more involved, 
although a significant and perhaps unnatural tuning of 
both the initial conditions and the potential seem neces- 
sary. 

Of course, care must be taken with the above argu- 
ments, as the energy exchange between the particles must 
be considered throughout the particles' orbits, and not 
just at the outset of their motion. Nevertheless, as we 
shall discuss further below, we have been able to ver- 
ify analytically that perturbations in nearly-harmonic 
potentials (to be defined more precisely below) with 
a;'((/)o) < are unstable, while those with a;'((/)o) > 
are stable. We will then discuss numerical results that 
support these conclusions. 



IV. QUANTITATIVE ANALYSIS OF 
DYNAMICAL INSTABILITY 

We now return to Eq. pTjl to discuss the quantitative 
behavior of perturbations. In a harmonic potential, V" — 
up' is constant, and (50 oscillates with fixed amplitude for 
all fc^. In the most general potential, V"{(p) oscillates — 
although not necessarily sinusoidally — with an oscillation 
frequency 2w, and so the perturbation equation, Eq. 
is that of a harmonic oscillator with an oscillating mass. 

The solutions to Eq. (fTT|) , as well as the issue of their 
stability, is the subject of Floquet theory. -'^ There is no 
simple stability condition in the most general case, but 



If we replace time 4 by a position x, our equation becomes the 
Schrodinger equation for a particle in a periodic potential, the 
solutions of which are described by Bloch's theorem. 



we can obtain an analytic solution for potentials that 
are close to harmonic. To begin, consider the quar- 
tic potential ^(0) = {\l2)n-?(\?- + (A/4)(/)^, working to 
linear order in A in the limit A 0. In this case, 
y'((/)) — TO^-|-3A0^. If the oscillation amplitude is 00, the 
oscillation frequency is uP' = m^[l 4- (3/4)A0o/m^] and 
the homogeneous oscillation is nearly sinusoidal: 0(t) — 
(j>ocosujt. With the trigonometric identity, cos'^ a; — 
(1/2)(1 -I- cos 2a;), the perturbation equation is, to linear 
order in A, then 

Scjj +{[k^ + 771^ + (3/2)A0^] + (3/2)A0§ cos2Lut} Scj) = 0. 

(12) 

Defining z = wt, this is identified as the Mathieu equa- 
tion, 

d? 

— (50-f [a-2gcos2z](50 = O, (13) 

with 

\ m'^ J \ 4 TO^ / 2 4 m'^ 

(14) 

If we consider only long- wavelength (i.e., k <^ m) fluc- 
tuations and the limit A(/>o ^ '^^j then a ~ 1 and 
q <^ \. In this regime, 5(j) oscillates rapidly with fre- 
quency Lot with an amplitude that varies as e^^* with 
~ -^q^ — (1 — af'LO. Thus, the condition for instabil- 
ity is 1 — |q| < a < 1 + \q\ (as derived in Section IVll 
below.) 

Thus, if A > 0, then the solutions are stable for all 
k^ > 0. For A < 0, the solutions are stable only if 
k^ > k^j — — (3/2)A(/)§. Recalling that the oscillation 
frequency is lu'^ ~ + 3A0§ , we see that there is insta- 
bility, on sufficiently large scales fc~^, if w'(0o) < 0, while 
the perturbations are stable for all k if a;'(0o) > 0. 

We can also generalize to other potentials that are 
close to harmonic, by which we mean that the time de- 
pendence of y"[0(i)] can be approximated as V"{t) « 
Vg" + V^2"cos2wt with constant Vq and V^a". The Jeans 
wavenumber kj is then given by 

k'j = l\V^'\+u;'-V^'. (15) 

It can be checked numerically that for the power-law po- 
tentials V{4>) oc l^l" with n ~ 2, these relations imply 
stability for n > 2 and instability for n < 2, as our 
heuristic arguments suggest. 

Of course, the types of potentials that can drive cos- 
mic acceleration are far from harmonic (for the power-law 
potential, an index n < 1 is required), and so the anal- 
ysis presented above will no longer be valid. However, 
we expect no qualitative difference, and to confirm this 
we investigate the stability of oscillations in more com- 
plicated potentials numerically in the next Section. 



V(^) 



V(0) 




FIG. 2: Two examples of a potential that can produce accelerated expansion and have stability on small scales — but only with 
oscillation amplitudes near those indicated by the dots. 



V. NUMERICAL RESULTS ON DYNAMICAL 
INSTABILITY 

We will numerically determine the stability of oscil- 
lations in the presence of three different representative 
classes of potential in this section: power-law potentials 
with arbitrary index, temporarily stable potentials ex- 
hibiting cosmic acceleration (as in the left cell of Fig. [2]), 
and stable potentials exhibiting temporary cosmic accel- 
eration (as in the right cell of Fig. [5]). In each case, we 
write the potential in the form V{(l)) — fj,'^v{(p/(j>o), and 
define the dimensionless variables. 



y 



where A is an appropriately defined constant character- 
izing the amplitude of the perturbation. The equations 
of motion then become 



= 







dT^ 



dv 
dx ' 



dx'^ 



y- 



(17) 



To determine the regions of stability, we can form the 
fundamental solution matrix (for a more detailed discus- 
sion of determining the stability of Hill's equation, see, 
e.g., Ref. M) 



C 



2/1 (T/2) 2/2(7/2) 
dyi/dT{T/2) dy2/dT{T/2) 



(18) 



where 2/1 is the solution generated from the ini- 
tial conditions {y{0) = l,dy/dT{0) — 0}, and 2/2 is 
the solution generated from the initial conditions 
{2/(0) = 0,dy/dT{0) = 1}. The stability of a solution 
can be determined by the eigenvalues of C (noting that 
det{C) = 1) 



Ah 



tr{C) ± VMC? 



(19) 



Instability results when either eigenvalue has modulus 
greater than unity. Since A+A_ = 1 (from the determi- 
nant condition), if the roots are both real, then solutions 



are unstable. This occurs when |i?'(C)| > 2. It also fol- 
lows that on the boundary between stability and insta- 
bility, yi{t) is a function of period T/2 or T (either the 
same period as v" or x). Because the system of equa- 
tions we are studying are invariant under time shifts by 
a period of oscillation, from 2/1 (0) — dy2/dT{0) — 1, we 
wiU have 2/1 (T/2) = dy2/dT{T/2). The stability criterion 
can in this case be phrased as \yi{T/2)\ < 1, and pro- 
vides an easy way to numerically determine stability or 
instability. 

We will first consider power-law potentials 



vm) = x\^\^ 



(20) 



The singularity at the origin can be regulated by intro- 
ducing a numerically small smoothing parameter c, such 
that 1^(101) = A(c-H^2)«/2_ i^g^yg checked that the nu- 
merics are insensitive to the value of this parameter. Nu- 
merically integrating Eq. ([T7)) . and solving for the roots 
of |2/i(T/2)| — 1, we can find the boundaries of stability. 
Plotted in Fig. [3] are the seven lowest bands of instability 
(the shaded regions of the Figure) for power-law poten- 
tials with n < 3. It can be seen that there is a band 
of instability encompassing k — for n < 2 as expected 
from the general arguments given above. In addition, 
as n decreases, the higher-order bands of instability mi- 
grate towards lower wavenumber. For power-law poten- 
tials with an equation-of-state parameter w ~ — 1, the 
power-law index will be rather small, and the large num- 
ber of higher-order bands of instability at small k might 
become important in analyzing the overall stability of 
such a model. 

Moving to potentials of the form shown in the left cell 
of Fig. [21 which we can write as. 



Cl 



C2X" 



(21) 



we chose as an example ci = C2 = 0.005 and n ~ 10. 
Again, numerically integrating Eq. (fTT]) , we plot the first 
six bands of instability in Fig.|4l Also over-plotted are the 
values of the equation-of-state parameter corresponding 
to various amplitudes of the background oscillations. It 
can be seen that there are indeed regions where w < — 1/3 
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10 15 20 



FIG. 3: The first seven bands of instability for oscillations in 
the presence of a power-law potential. 




FIG. 4: The first six bands of instability for oscillations in the 
presence of the potential shown in the left cell of Fig. [2] The 
solid line represents an amplitude corresponding tow = —0.4, 
dashed line w — —0.6 and dot-dashed line w = —0.8. In this 
example, w remains less than one, but as the amplitude of 
oscillations decays a large-scale band of instability develops. 

and large-scale stability exists. However, as the universe 
evolves and the amplitude of oscillation decays, more 
and more bands of instability will pile up at small k. 
When the amplitude of oscillation dips below xo ~ 1 in 
this model, the solution will become violently unstable 
at both small and large k, with only very small bands 
of stability. This will most likely present a severe chal- 
lenge to using such models as a candidate for the dark 
energy. We have also checked numerically that the large- 
scale (encompassing fc = 0) band of instability appears 
when the amplitude of oscillation is such that w' changes 
sign, supporting, in a more complicated potential, our 
heuristic arguments for stability. 

Finally, we will treat the case where we expect the 
potential to exhibit large-scale stability, but only tempo- 
rary cosmic acceleration, as in the right cell of Fig. H 
The potential is given by 

V{cl,)^fi^x^ + c')\ (22) 

where we choose as an example c = 1. When the am- 
plitude is very close to the height of the barrier between 



FIG. 5: The instability band for oscillations in the double- well 
potential shown in the right cell of Fig. [2] There is just one 
band of instability until k becomes very large. Also plotted 
are lines that indicate the amplitudes at which w = 0.4, w = 
0.2, and w = (from top to bottom). Note that w ^ — 1 
as a; ^ \/2 ~ 1.41, and so this model features large-scale 
stability with — 1 < w < 0. However, this period of negative 
pressure is short-lived as the Universe expands. 

the two minima of this potential (at Xg — \/2), wc ex- 
pect the equation-of-state parameter to be ~ — 1 ^. As 
shown in Fig.O this potential yields large-scale stability, 
with the first order band of instability approaching fc = 
as the amplitude of oscillations reaches xq ~ ^/2. Also 
shown are the values of w for a number of amplitudes. 

In all of our examples, we have verified numerically 
that the sign of uj' determines the existence of a band of 
instability about A; = 0, as expected from the heuristic 
arguments of Sec. IIII Bl We have also illustrated the 
existence of higher-order bands of instability that will 
infiuence the development of small-scale inhomogeneities 
in the oscillating field. In all examples, the higher bands 
of instability creep in towards fc = Oasu'^— 1. 



VI. GRAVITATIONAL INSTABILITY 

A. General analysis 

We now include gravity in the analysis. As alluded to 
above, perturbations in a perfect fluid with sound speed 
Cs are stabilized by pressure gradients for wavenumbers 
k > kj = {A-KGp/ciy^'^, but care must be taken in ap- 
plying this result to coherent scalar fields. For exam- 
ple, oscillations in a harmonic potential have a sound 
speed = d{wp)/dp = 0, implying that perturba- 
tions on all scales will suffer a gravitational instability. 
However, as shown in Ref. [l^l, perturbations on scales 



^ The field will loiter near the local maximum, extending the pe- 
riod of oscillations, and eventually spoiling our approximation 
that the period be shorter than a Hubble time; we will ignore 
this shortcoming for the moment. 
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k > 2v7rGTO^0o are stabilized. We therefore extend the 
analysis of Ref. [13] to see how their result for the scalar- 
field stability is altered by the presence of a small an- 
harmonic term in the potential. This is in fact general, 
since, as we will see, strong dynamical instabilities will 
completely swamp the gravitational instability, implying 
that the inclusion of gravitational effects is only relevant 
for potentials that are close to harmonic (as we have de- 
fined above). 

To simplify the analysis, we restrict our attention to 
small-wavelength Fourier modes, fc 3> iJ, where H = 
(87rGp/3)i/2 and p = (1/2)02 ^ y^^y j.^^^ j^y^^g ^3 

to neglect the expansion and work with a perturbed 
Minkowski metric (in the conformal Newtonian gauge), 



(1 + 2'^)dt^ + (1 - 2^)dx^. 



(23) 



where ^' is the Newtonian potential. Perturbations Scj) 
in the scalar field may now induce perturbations '5 to 
the metric, which then affect the scalar-field equation of 
motion, which is now, 



'+[v"m)) 



^ m{S(t>)-2^V' {(l}{t)) . (24) 



The potential is determined by the Einstein equation (the 
Poisson equation). 



V (0(i)) (5(/. - 
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(25) 



We now focus on a given Fourier mode of wavenumber k, 
in which case we replace — > —k^ in these equations, 
and we neglect the last term on the right-hand side of 
Eq. ((25|) . as it is negligible ior k ^ H. 

Although the equations have become more complicated 
with gravity, they are once once again linear differential 
equations with periodic coefficients, and so the solutions 
are formally described by Floquet theory. 



B. A quartic correction to a harmonic potential 

We now consider the potential V{4>) = (1/2)™^!/)^ -t- 
{X/4)(f)'^ and suppose that the anharmonic term is small, 
X4)l < Eqs. (I24|) and ^ can then be combined 

into a single second-order differential equation. 



S4> + 7717 sin 2ujt5(j) + 



fc^ + + ^A(/)o - ^ni^ + f ^'^'/'o 7"^^ 1 cos2wi 



(26) 



where 7 = Gcft^/k"^. At the lowest- A: stability boundary, 
the solutions are periodic with frequency w, and can thus 
be written 5(t){t) — a cos ujt + bsinujt. We then plug this 
into Eq. ((26)) and demand that the coefficients of the 
terms that vary as cos ujt and as s'm ujt vanish. In doing 
so, we work to linear order in A and 7 and neglect terms of 
0(A7). We then find that the boundary between stability 
and instability occurs for wavenumbers fci and ^2 given 
by 

kf^-hcPl+jm^, kl = 0. (27) 

A bit more algebra (replacing a and b by slowly varying 
functions a{t) — aoexp(rii) and b{t) = baexp{flt), with 

<C w) shows that the instability occurs for values of k 
between ki and ^2. 

If A < 0, then kf > 0. In this case, the field is dynam- 
ically unstable, even without gravity (as noted before), 
and gravity only serves to increase the instability. If, 
however, A > 0, gravity induces an instability (noting 
that 7 cx fc^2), for k^ < kj, with 



k'j^-l\c^l + ]J(jX^l^ +87rGm402. (28) 

Note that this recovers our earlier result, without gravity, 
in the limit that G ^ 0, and it recovers the result of 



Ref. [23,[25^ for A = 0. 

The result for the Jeans scale resembles that for 
spintessence (cf., Eq. (10) in Ref. Q) with the replace- 
ment {l/2){V'/\(j)\ - V") = -A|0g| for spintessence with 
— (3/4)A0g for the oscillating field, and replacing the G 
in the spintessence result by G/2. The former replace- 
ment occurs because the sound speed cj. = d{wp)/dp = 
(l/2)A(0Q/m2) for spintessence differs from that, = 
(3/8)A(0q/to2), for the oscillating field. The latter re- 
placement (i.e., G G/2) occurs because the com- 
plex field is equivalent to two scalar fields. Note that 
the result, Eq. ([28]) , differs from the spintessence re- 
sult also because the spintessence Jeans scale of Eq. (10) 
in Ref. [1] places no restrictions on the potential V{R), 
while Eq. here is valid only in the limit A(/)q <C m? ■ 

In the hmit that Gm'^^g < \^(t)l, kj ~ 167rGTO*/3A, 
which is equal to 47rG/c^, the Jeans scale for a perfect 
fluid. However, if Gm^(/)§ ^ A^t/ig, then the perfect- 
fluid description breaks down, the scalar-field dynamics 
become important, and the Jeans length differs consid- 
erably from the perfect-ffuid result, c^/cj = AirGp. 

C. General result for nearly harmonic potentials 

The result, Eq. ([^5)1 . can be generalized to other poten- 
tials that are close to harmonic by replacing —{3/2)X4)q 
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in Eq. ^ by (1/2) IK,"! + - V^', 




V" 



+ 87rGm402 



.(29) 



range 3000 ^ 300, resulting in a small-scale cutoff of 
order 10~^^ in the primordial power spectrum, much 
smaller than that for WIMPs. Note that there may also 
be dynamical instabilities on even smaller scales (related 
to the higher-order bands of instability, as in Sec.|V|; see, 
e.g., the very elegant work of Ref. [33l |. 



D. Application: Axion dark matter 

We now calculate the small-scale cutoff in the cold- 
dark-matter power spectrum under the assumption that 
the dark rnatter is composed of axions with masses ma ^ 
10~^ eV [29]. This cutoff will determine the masses of 
the first dark-matter halos to undergo collapse in the Uni- 
verse, and it determines the size of the smallest clumps in 
the Milky Way halo [s^, [3l| . If weakly- interacting mas- 
sive particles (WIMPs) make up the dark matter, then 
the primordial power spectrum is suppressed for cosmo- 
logical mass scales smaller than ~ 10~^ — 10^ (the 
precise value determined by the precise WIMP model) 
|33 | by kinetic decoupling of WIMPs. We will now cal- 
culate the analogous small-scale cutoff if axions make up 
the dark matter. 

Cosmological axions in the ~ 10~^ eV mass regime 
are produced by the misalignment mechanism near the 
time of the QCD phase transition, and they may thus be 
described by a coherently oscillating scalar field. The 
axion potential has the form V{(t)) = Vo[cos(0//) — 
1]. Today, the axion field oscillates near the mini- 
mum of this potential where it can be approximated by 
V{<j)) « (l/2)m202 _^ (A/4)04^ with m2 = Vq/P and 
A = -(l/6)(mV/^)- Here, / ~ A^/m is the Peccei- 
Quinn scale, where A ~ 100 MeV is the QCD scale. If 
axions make up the dark matter, then the axion energy 
density today is pq = {l/2)m?(j)l — flcPc where fie ~ 0.2 
is the cold-dark-matter density and pc is the critical den- 
sity, and this relation can be used to fix 0o. Using the 
scaling oc (1 -I- z)^, we find also that (Ai/)^)^ ^ Gm^cj)^ 
not only today, but at all redshifts z < 3000 during mat- 
ter domination. The axion-axion interactions implied by 
the A(/)^ correction to the quadratic potential therefore 
have little effect on the Jeans scale, which is well approx- 
imated by the earher result of Refs. [2^[25t. 

The physical Jeans scale is given by kj — 2^/2tiG7ii?4>q, 
and the Jeans mass then turns out to be Mj = 
(47r/3)(7r/fcj)^p„j, where pm is the current matter den- 
sity. Numerically, 

Mj ~ 1.8 X 10~" (m„/10-5 ey)-3/2(l + z)-^/'^M<s. 

(30) 

The collapse redshift for such low-mass halos depends 
on the primordial spectral index, but is generally in the 



VII. CONCLUSIONS 

We have considered the dynamical and gravitational 
amplification of inhomogeneities in oscillating-field mat- 
ter. We provided a simple physical picture for the 
origin of instabilities in oscillating-field matter and in 
spintessence. In this picture, the condition for instabil- 
ity is seen as a condition on the amplitude dependence 
of the frequency, a condition that can be shown to be 
equivalent to positivity of the sound speed. This argu- 
ment was verified analytically for nearly harmonic po- 
tentials and numerically for more general potentials. We 
then included gravity in the analysis, generalizing earlier 
results on the gravitational instability of harmonic po- 
tentials to potentials that differ slightly from harmonic. 
We used this result to evaluate the small-scale cutoff in 
the matter power spectrum if axions make up the dark 
matter. We leave the implications of this cutoff (roughly 
10^^^ M©) to future work. 

Our results indicate that potentials that give rise to ac- 
celerated expansion generically suffer dynamical instabil- 
ities to the growth of large-scale inhomogeneities. These 
instabilities should render oscillating fields unsuitable to 
account for dark energy in the Universe today or for driv- 
ing infiation in the early Universe, as both scenarios re- 
quire the cosmological density to remain homogeneous for 
extended periods of time. It is true that there may be 
potentials that drive acceleration and are stable, but sta- 
bility and/or acceleration will be only temporary. This 
loophole is thus unlikely to alter our conclusions. 

We were able to make progress analytically for both the 
dynamical and gravitational instability only for nearly 
harmonic potentials. It would be interesting to see 
whether analytic results for the growth of perturbations 
can be extended to more general potentials. We leave the 
investigation of this question for future work. 
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